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1. Introduction 



Data from heavy-ion collisions from RHIC and LHC have revolutionized our per- 
ception of strong coupling physics in QCD, and revealed the characteristics of the 
deconfined phase. They have also become the testing ground of novel techniques 
emerging in string theory that attempt to control strong coupling phenomena using 
a gravitational description. In this description, a heavy-ion collision is described as 
the process of black hole formation and decay, albeit in a five- dimensional theory of 
gravity including also other fields, notable a scalar, the dilaton, fl]]. 

In this direction, several attempts have been made to analyse the scattering of 
high-energy sources, using shock waves in AdS, [H, ||, |], ||, ^, ||, |10], jlT], [12], |13 

T|, |i5|, |i6|,|i3,|i8i||, igi. 



A direct outcome of this approach, is the estimate of lower bounds of the final 
multiplicity by using Penrose's idea of trapped surfaces. 

In this paper, we analyse the formation of trapped surfaces in head-on collisions 
of shock waves in gravitational theories with more complicated bulk dynamics, and 
different types of vacuum solutions. An example are Einstein-dilaton theories with 
a scalar potential. Such theories have been argued, p6|, to describe holographic 
physics that is closer to QCD than the AdS theory. A phenomenological theory. 
Improved Holographic QCD (IHQCD) has been constructed, |2^, |2^, that agrees 
well with both zero temperature and finite temperature YM data. 

The idea is to explore how different aspects of the dynamics affect the trapped 
surface that forms during the collision of shock waves. The shock waves used are 
gravity (spin-2) shock waves. The different factors are the following: 

• Different bulk geometries. There are several different possible geometries that 
that have been classified in [26||. They are characterized by their IR and UV 



(near-boundary) behavior. 

In the UV the typical behavior is asymptotically AdS. There is however inter- 
est in different non-AdS asymptotics, as they can capture the physics of the 
collision when the trapped surface forms mostly in the IR part of the geometry. 

The IR behavior can be split into three large classes: (a) Confining geometries, 
(b) Non-confining geometries, (c) Unacceptable geometries (that violate the 
Gubser bound^). AdS^ belongs to class (b). There is another special geometry 



^Thcsc works refer to the AdS^ geometry whose dual theory is the J\f= 4 sYM and not QCD. 
References |2^, review the similarities and differences of the two gauge theories. The works 
pl| , p2| , ^ discuss the recent development in the field of applications of AdS/CFT in QCD. 
Particular emphasis in heavy-ion collisions and the Quark-Gluon Plasma is given in [ p^ p5[ . 

^This is a criterion introduced by Gubser in in order to test the acceptability of solutions 
with an IR singularity. It states that a naked singularity is acceptable ( "good" ) if it can be covered 
by an infinitesimal horizon. The implications of this constraint for Einstein dilaton gravity were 
analyzed in p6[. 
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that lies at the boundary of (a) and (b), the hnear dilaton geometry. 



Different profiles of the transverse distributions of energy. So far two types have 
been analyzed, uniform distributions or power-like (aka GPY, [^]) profiles. We 
will add one more class namely exponential profiles, that are well localized in 
the transverse plane. Such energy profiles are closer to what the targets are 
in heavy-ion collisions. Except the uniform distributions, the other profiles are 
characterized by a length scale that controls the size of the energy distribution. 
Typically GPY profiles are fuzzier as they stretch to larger distances, while 
exponential profiles are sharply localized. 



• The option of cutting off the UV part of the bulk geometry, ||3T|]. This is mo- 
tivated by the fact that we are striving to emulate the QCD behavior that is 
perturbative in the UV. In a perturbative regime, by definition, multiplicities 
are small. It is in this regime that the geometry is expected to be unreli- 
able in a holographic description. The simplest way to implement this is to 
stop the geometry at a point deemed to be the transition to the perturbative 
regime, and ignore the contributions above that point. This is not a controlled 
approximation, but it is expected to give useful hints. ^ 

Our approach is to solve first the equations for the shock waves, in the presence 
of different types of metrics, with sources in the bulk. Such solutions determine also 
the transverse distributions. We then find the associated trapped surfaces, and the 
estimated lower bound on the generated entropy'^. We will use the word entropy 
liberally in this paper to mean the area of the trapped surface. Strictly speaking 
this is a lower bound on the entropy, but it will always be obvious which entropy we 
are referring to. Moreover, entropy can be converted to multiplicity, and it is in this 
sense that we will use this two terms interchangeably. 

Two characteristic examples are analyzed quantitatively and eventually com- 
pared to data. 

1. The first is AdSs, with a transverse distribution having a localized exponential 
profile and a cutoff in the UV at r = l/Qs with Qg the saturation scale (section 
p.2.2| ). We will name this setup "AdS-Qs" for future reference. 



2. The second is a simplified metric emulating the IHQCD solution of [26 



im- 



plementing confinement and asymptotic linear glueball trajectories, while it is 



^In the bottom-up models of the gauge couphng constant, epitomized by the exponential 
of the dilaton is becoming small in the UV region. Despite this, the interactions of spin-2 matter 
are large in the UV. A potential way out is to advocate an asymptotically AdS metric in the string 
frame, but this cannot be done with an action with two derivatives, |65|. 



^An interesting approach for computing the entanglement entropy |32 and probing the scale- 
dependence of thermalization is provided by reference p3|. 



3 



asymptotically AdSs in the UV. In this case we have again a localized expo- 
nential profile in the transverse plane and no UV cutoff, (section |8.2.3| ). We 
will name this setup, with a slight abuse of language, "IHQCD" . 

What we find is as follows: 

• Both, our analytical calculations and our numerical analysis show that most of 
the entropy comes from the UV part of the geometry provided that the geom- 
etry reduces to the AdSs space near the boundary. This is not in contradiction 
with the naive expectation that in a strongly coupled theory the multiplic- 
ity comes from the high-energy part of the phase space available during the 
collision. 

• For uniform transverse distributions, the AdSs geometry, produces the least 
entropy among the geometries of section |^ . For non-confining geometries the 
entropy scales as for large s, with | < a < | with a = | for AdSs. For scaling 
confining geometries | < a < 1. The case a = | corrected by logs corresponds 
to the geometries related to IHQCD. We have assumed above that the other 
geometries persist up to the boundary, or equivalently the trapped surface 
forms in the regime in which they are valid. These results are summarized 
in table 0. However, as most realistic geometries are asymptotically AdS, the 
thrust of this result may be important only at intermediate energies, where IR 
geometries may take over. 

• For geometries that have a mass gap, discrete spectrum and confinement, the 
allowed transverse distributions has a spectrum of scales that is in one to one 
correspondence with the discrete spectrum of 2++ glueballs. The associated 
entropy production, independently of transverse distribution, is less than AdSs. 
It is assumed that all such confining geometries are asymptotically AdSs. 

• An important puzzle of our analysis in the previous item is that the scattering of 
a distribution with transverse size associated with the lowest lying 2^^ glueball 
does not seem to lead to a trapped surface. 

• There are geometries with a UV energy-independent cut-off that lead to an 
asymptotic ~ log^(s) behavior for the entropy (see table 

• A general trend is that at equal total energy, the collision of distributions 
that have a larger transverse size leads to a larger entropy production. This 
implies that more dilute energy distributions produce more entropy at fixed 
total energy. In particular a uniform transverse energy distribution produces 
(at equal total energy) more entropy that one with power-law or exponential 
transverse distribution. Similarly a power-law transverse distribution at the 
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sane total energy produces more entropy compared to an exponential transverse 
distribution with the same length scale. 

• Comparing an AdS geometry, with an asymptotically AdS geometry that is 
confining in the IR we find that for all different types of transverse energy 
distributions that we have examined (uniform, power-law, or exponential), the 
AdS geometry generates substantially more entropy than the confining geom- 
etry, at the same (and large ) total energy, and transverse scale. 





b{r) 


L 

r 


-exp[-;^] 


Transverse profile 








Uniform 




cAdS 


nlHQCD 
unif 


GYP 




cAdS 
'-'GPY 


Not studied 


Exponential 




cAdS 
^exp 


qlHQCD 



Table 1: The several cases analysed and compared. Two geometries have been compared: 
A non confining geometry (the AdS^) and a confining [IHQCD) one. None has a cutoff. 
The transverse profiles correspond to the cases of being constant, GPY (falling-off as a 
power law, [pl|) and falling-off exponentially. 



More to the point, table |l] defines the different contexts studied in this paper. 
We find the following inequalities between the various trapped surface areas 
(keeping total energy and transverse size fixed) 

oAdS ^ cAdS ^ cAdS qIHQCD qIHQCD /i i\ 

'-'unif ^ '-'GPY ^ '-'exp J '-'unif ^ '-'exp {^■^) 

qAdS ^ qIHQCD qAdS ^ qIHQCD qAdS ^ qIHQCD /i 

'-'unif ^ '-'unif ' '-'GPY ^ '-'exp ; '-'exp ^ '-'exp K^-^) 

In confining backgrounds, the entropy production increases as the confinement 
scale Aqcd decreases, provided that the total energy and the transverse size 
are kept fixed. 

In the AdS-Qs setup the multiplicities grow with the atomic number almost 
linearly; in particular as ~ 



• The lower bound on multiplicities for AdS-Q^ is given in equation ( |8.4| ). For 
the IHQCD setup it is given in for the (|8.6|) . Both formulae once fit to RHIC 
data, make the same (correct) prediction 2.76 TeV PbPb LHC data, [^. A 
single parameter is used in these fits, namely the overall constant coefficient of 
the leading (large) s-dependence. 

The implications of the conclusions above do not a priori apply verbatim to QCD 
because: 
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1. We are discussing glueballs only, and not nuclei. 



2. The holographic backgrounds used are approximations to various gauge dynam- 
ics. Only AdS is an exact background for sYM, but all backgrounds, including 
AdS used for QCD, are approximations. 

3. The final entropy produced may be and is usually larger than the trapped 
surface area. 

The first issue may not be of prime importance at high energy. The reason is 
that the energy released at mid-rapidity in heavy-ion collisions is expected to be 
mostly gluons. Quarks will provide corrections to this but they are not expected to 
change this picture drastically. 

The second issue is of importance. Although by now bottom-up models of YM 
can provide reliable calculations matching lattice calculations, their structure in the 
UV is more shaky. It is not a priori clear where the geometric description breaks 
down, and such transitions may be at different places for different observables. 

The third issue is also important for quantitative predictions. There are very 
few cases where the final entropy has been calculated numerically in the collision of 
shock waves, JIT], . In [0 in particular it was shown that the released entropy at 



the end of the collision process is 60% larger that the bound found from the trapped 
surface calculation. Moreover conformal invariance in AdS makes this percentage to 
be independent of the collision energy. 

This suggests that for collisions and bulk geometries where the majority of the 
trapped surface area come for the UV, AdS part of the geometry, we should still 
expect that the relative factor relating the final entropy release to the area of the 
trapped surface to be almost energy-independent. 

The analysis in this paper is providing important and potentially general clues 
on multiplicity generation in high energy collisions using holography. It must be 
however be backed-up by a more reliable calculation of the gravitational evolution, 
but we leave this for future work. Another interesting issue is whether there are 
differences between the high-energy scattering of 2^"*" glueballs studied in this paper 
and 0^"*" glueballs. Although we do not have reasons to expect major differences, an 
analysis in this case should be done. 

We organize this write-up as follows: We begin in section ^ by showing the way 
one may build more realistic geometries that come closer to QCD. This is achieved 
by considering scalar gravity. In section 0, we state the equations that compute 
the entropy production of the collided shock-waves and hence of the gauge matter 
indicating that the calculation reduces to a boundary valued problem. In section ^ we 
present the subclass of geometries we investigate. Sections |^ and | refer to various 
geometries with uniform or non-uniform transverse dependence respectively. The 
entropy of the trapped surface is computed for each case. Section ^ deals with the 
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simplest way of removing the weak coupling entropy production: At higher energy 
scales, where coupling is weak, the contribution to the entropy production should 
be less important. In the geometrical language of AdS/CFT this would imply that 
one should modify the way the fifth coordinate of AdS/CFT at the UV is treated. 
In section |^, we present our results and make our predictions. Finally, in section 
1^, we comment on various aspects of our investigations and conclusions. We omit 
intermediate details of our calculations for the three appendices at the end. In 
particular, appendix ^ serves as a practical introduction to the theory of trapped 
surfaces. Appendix |B| shows how one may localize the (five-dimensional) bulk sources 
and appendix proves a useful equation that is needed for subsection |6l^ . 



2. Einstein-Dilaton gravity 

We start from the action: 

S5 = -M' j (fx^ 



(2.1) 



where $s is the scalar field dual to the YM coupling constant. We first find shock 
wave solutions in this theory of the form 

ds^ = h{rY [dr"^ + dx'-dx'- - 2dx^dx~ + (f){r, x^, x^)5(a;+)(dx+)^] , = ^sir, x+) 

(2.2) 

with the asymptotically AdS boundary at r = 0. Compatibility of these equations 
implies that = 0. Eliminating $5 using the equations of motion, the equation 
for 6 is 



b' 



Vi + 3-dr + dl 
b 



= -24J++, Vl=didi, b' = drb{r) 4 = 87rG5 (2.3) 
where we have introduced a stress-tensor 7++. 



3. Trapped surfaces 

Trapped surfaces are created when two shocks like the one of ( |2.2|) ^ which moves 
along x~ collide^. In terms of metrics before the collision, one then has 



ds'^ = b{r f dr^ + dx'dx' -2dx^dx + (f)i{r, x^ , x^)5{x^){dx^y 



(3.1) 



+ (f)2{r, x^ , x'^)S{x ){dx 



x± < 0. 



^We call it 0i to distinguish it from the second that moves along x'^ that we call 4>2- In our case 
the shocks will be taken identical and hence the subscripts will be soon dropped. 

more complete set of notes on the theory of trapped surfaces may be found in appendix 
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Associated with the shock-wave (pi in (|2.2|) , we parametrize (half of the) trapped 
surface 5*1 by 

x+ = X- + ^i:i{x\x^,r) = (3.2) 

where i/^i remains to be determined. It is also useful to rescale the functions 0i and 
ipi by defining 

$1 = b{r)(j)i ^1 = b{r)^i. (3.3) 
\l'i satisfies the following differential equation 

40 = A. MM (3.4) 

where O^^g^ is defined with respect to the metric 

ds^ = h{rf [dxl + dr^) , dxl = {dx^f + {dx^f. (3.5) 



We point out that once ( ^.3| ) is solved, then (DAdSs ~ ■^ip{r)))^i provides a source 
term for (ClAdSs ~ ^(^(^)))^i- The missing ingredient is the boundary conditions 
that are given by 



^1 



= E [V^^i^^^i] =8 (3.6) 

c h[r)^ c 



i=l,2,r 



for some curve C which defines the boundary of the trapped surface and where both, 
Si and (the associated surface to 02,) 'S'2 end. The entropy is then bounded below by 
the area of the surface obtained by adjoining the two pieces of the trapped surface 
associated with each of the shocks as 



S>Strap = '2x^ / ^Jdet\gAdsMrd?x^ = ^ I Krfxl{r) (3.7) 

4tJ5 Jc ^'-^5 Jrc^ 

where the (generalized) curve C defines the boundary of the trapped surface Si and 
S'2 which are identical; thus the overall factor of 2. The integral with respect to the 
transverse coordinates gives x\{r) when considering a head-on collision. Typically, 
f'ci, ^± carry the information of the shock ''. We conclude by pointing out 

that the un-scaled version of the trapped surface (see (|3.3|) ) for a head-on collision 
is given by 

Vl + 3^9, + 9,2^ (0-V) = 0, ^ = 0|^, (drijf + {d,^^)\ = 8 (3.8) 

while the entropy is still given by ( p.7| ). In the absence of transverse dependence one 
ignores x\ from ( p.7| ) which measures entropy /transverse area in this case. 

^We drop the subscripts 1, 2 from 0's and V^'s from now on. 
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4. Shock geometries 



We have analyzed a number of bulk geometries by examining a class of different scale 
factors, that classify the non-conformal behavior of Eistein-dilaton gravity models. 
In terms of their behavior in the far IR (r — > oo, or r — )■ tq), we have the following 

cases 



26, m 



1. b r"", with a < —1. The AdS case corresponds to a = 1. This corresponds 
to quasiconformal geometries, with no confinement, continuous spectrum and 
a mass gap, with potential asymptotics as $s — t- oo, V ~ e*^*% Q < |. 



2. Confining backgrounds that are scale invariant in the IR, |56|, with b{r) ~ 
(ro — r)", a > |. In this case tq is finite and signals the position of an IR 
singularity that satisfies the Gubser bound for a > |. They have a discrete 
spectrum of glueballs and a mass gap. The potential asymptotics as $s — )■ oo 



are V ~ Q > 



3. Confining backgrounds with b{r) ~ e (^'')'', a > 0. They have a discrete 
spectrum and a mass gap. The potential asymptotics as $s — ?■ oo are V ~ 

4 A S—1 
6 3*= $ " . 



4. Confining backgrounds with b(r) ~ e V'— '■oy , a > and Tq as in the second 
point above. They have a discrete spectrum and a mass gap. The potential 
asymptotics as $s — ^ oo are V ^ es^" ■ 

5. Shocks with uniform transverse space dependence 

In this case the shock (f){x~^,r) can be determined and the trapped surface ends at 
r = Th with Th determined from^ 

/rlr F 
^ , b'{rH) = -^ (5.1) 

with E s^, and s the center of mass energy of the collision. The area of the 
trapped surface is 

Atrap^ I " b'dr (5.2) 



We may therefore estimate the energy dependence of the trapped area for different 
bulk geometries. 

• For non-confining scaling theories b ~ r^", with a > 1 we obtain Atrap ~ . 
The AdS case corresponds to a = 1. This agrees with previous estimates for 
the AdS case. 



The delta function in (|5.lD may be in principle replaced by any function of . 
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• Confining backgrounds with 6(r) ~ (rg — r)", a > |. In this case we obtain 
Atrap ~ at high energy. The exponent varies between | and 1. 

• Confining backgrounds with b{r) ~ e"'-'^''''". In this case Atrap ~ s^(logs)~^ 
at high energy. 

• Confining backgrounds with 6(r) ~ e V'-'-o/ . In this case Atrap ~ s 2 (logs) « 
at high energy. 

Taking into account that Strap ~ Atrap we conclude that in all cases the entropy 
production is larger than AdS at high enough energies. 



6. Non-uniform transverse dependence 

We assume the shock(s) have non-trivial transverse dependence and we solve the 
homogeneous equation ( |2.3| ) by separating variables. We obtain the following set of 
differential equations (assuming rotational symmetry on the transverse plane) 



(6.1; 

The first equation yields 



A(x±) = CiKoikx^) + C2loikx^) (6.2) 

while the solution to g{r) depends on the scale factor b{r). Strictly speaking, the 
solutions corresponding to Kq do not exactly solve the homogeneous differential 
equation for f^. They induce a delta function source and they satisfy 

(vl + 3^dr + d^^ Mkx^)gk{r) = -27iS^^\^^)gk{r). (6.3) 

It is evident that J++ from ( |2.3| ) consists from an appropriate linear combination of 
KQ{kx±)gk{r)'s. In the remaining of this section we analyze several 6(r)'s. 



6.1 Power-like transverse distributions 



We will study solutions to (|6.1| ) where the transverse distributions fall-off as a power 
in the transverse plane. Such solutions were considered for AdS in ||31|. In the 
remaining of this section we analyze several different types of scale factors b{r) and 
we derive the boundary that specifies the trapped surface as a function of the energy 
(see equations (|6l7| ) and ( |6.13| )). Then using (3.6) one may compute the entropy 
bound S; the final answers for S are summarized in section 8 (see tables 0, 0, 11). 
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6.1.1 Power-law b{r) 

We begin by studying the power-like scale factor 

b{r) = ^{r-ror. (6.4) 

When a > I and tq finite, this corresponds to a confining geometry, with IR singu- 
larity at r = ro, This can be resolved by uplifting to higher dimensions, [^. It 
can also include non-confining geometries, if ro = and a<— 1. a = — 1 corresponds 
to AdS. The boundary for a>l/3isatr— )-oo while for a < is at r —t- 0. 

Assuming a point-like 7++ = E6{r — r')6{x^){x^)6{x^) (r' is the r coordinate of 

the source) in (|2.3|) , the shock is given by 

i_ 

0=(— — ) 2(l+3«)$(g), $(g)^ „_ ± V ; 



(g(l + g))^ ' 4(r-ro)(r'-ro) 

(6.5) 

where P^^ is the associated Legendre polynomial and where the tilde on $ is to remind 
us that the rescaling of (j)^ is not by a factor of b{r) to the first power (see (|3.3| )) and 
will be omitted from now on. The presence of the combination q, natural in AdS 
space, [^] for more general power-like metrics, can be explained, as such metrics are 
conformally AdS, and can be lifted to AdS metrics in a higher dimensional space- 
time, . 

The equation for the trapped surface satisfying the first boundary condition is 
given by 

$(gc) Q!i+3,(l + 2g) 

"^-{qc) (g(i + g))4 

where $_ is the solution to the homogeneous differential equation satisfied by $ 
which is finite at g = 0. is the associated Legendre polynomial of the second kind 
while qc denotes the value of q at the boundary of the trapped surface. The second 
boundary condition in ( p.8| ) yields 



^(l+3a) 



|(l+a) 7-1 2/ / \ M'c ! > q) 



'''' qcil + qc)^-ia;qc)^ { (6.7) 



-f(«-i) ^ 



where the quantity on the right of the arrow denotes the high energy limit. We note 
the agreement with when a = —1, vq = and r' = L. The Strap is given by 

Sa>i ~ / (r-rof^ Air - r,){r' - r,) -3-— - (r - rf) \ dr, 

(6.8a) 



^And likewise, the rescaling of -0. 
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Sa<-i ~ vr / (r-ror \ 4(r - ro)(r' - ro) " " ^')') | ^r. 

(6.8b) 

The limits of integration are found from ( |6.7| ) when x± = 0. Performing the integra- 
tions and keeping the leading contribution in E taking into account the range of a 
for each case, we finally arrive at 

^„>i/3 ~ (6.9a) 

(3o+l) 

^„<_i ~ E^^. (6.9b) 

As a cross check we note that ( |6.9b|) for a = —1 reproduces the result of [0. We 
remark that the trapped surface for a > 1/3 exists because tq is finite. 

6.1.2 Exponential b{r) 

We will take the scale factor b(r) in this case to be 

b = e-'/^. (6.10) 



This corresponds to a marginal confining geometry, that of a linear dilaton, 
Rescaling (p = eai? $(u(r, x_l))^° and assuming a point-like J4.4. as in the power-low 
case, we obtain 

$(n)^^-^e .v^, u^'- (6.11) 

where r' is the r coordinate of the source. The trapped surface equation is 

= - *_(,.) = (6,12) 

while the rescaled boundary condition -^e^'"/'^u(^'(M))^|^ = 8 finally yields 



E ~ 4^e-=^'^^/2«v^sinh(^v^) (6.13) 



where the quantity on the right of the arrow denotes the high-energy limit. The 
trapped surface is found by working as in subsection |6.1.1| . In this case also, the 
trapped surface exists either when r is allowed to extend in the interval [—00, 00] or 
when a UV cut-off is placed (see discussion below ( |6^) ). 

Allowing for r to take negative values, so that the geometry is extended towards 
the boundary, the entropy is computed from ( |6.13|) numerically: One, solves this 
equation with respect to x\ = x5^(r, E) and integrates as in (|3.7|) . The two limits of 
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And likewise, rescaling ip. 
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integration are found numerically from ( |6.13 ) again by setting x_l = 0. The analysis 
shows that 



Strap-- s'-^'\0g'-"is). (6.14) 

The result indicates a larger growth of Strap than the uniform profile case (see third 
bullet of section |^). This, does not contradict our general conclusions as here the 
exponential metric is assumed to be valid up to the UV boundary. There is a way to 
derive a simple formula for Strap when a UV cut-off (see section 0) is placed. If for 
instance, we cut the surface at some r > then equation ( |6.13| ) implies that as E 



increases, uc increases logarithmically: log'^{E) ~ + x]_q. Hence, equation (|3.71 ) 
yields 

5'tr-ap ~ ~ log^(s). (6.15) 

This result looks similar with Froissart bound estimates for cross sections p, ^ . We 
do not know if there is a connection. 

6.2 Localized transverse distributions 

We built a single-mode shockwave (j)k out of the solutions of ( |6.1| ) demanding to be 
square integrable for all x±. The result is 

— ^ jgi{kr)Ko{kx±)5{x^), assuming gi{kr) ^ (6.16) 



where E is the energy carried by the shock. fi'i,2('"), the radial solutions, are dimen- 
sionless, because (pk has dimensions of length. The overall constant ensures that the 
gauge stress tensor, which according to the AdS dictionary is given by 

T++ = ^ lim (6.17) 

when is integrated in space, it yields the total energy E. The trapped surface is 
defined by 

Vl + 3^dr + d^^ {(Pk - ^k) i^k = 0\^ {dr^kf + {d,^tPkf\c = 8 (6.18) 
and yields 

, , Ek' jClMk'r) + CjMk'r)) Ipjk'x^) ^ , 
^k = Mr,x±)- — — — ^ — — — rvYTT^ zMrc,x±c)- (6.19) 

for some coefficients Ck to be determined. The Kqs do not participate as they 
induce source terms^^. If we assume that at small r, gi — t- n(r) with n(r) a 
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In appendix H it is shown how these sources may be locahzed. 
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regular function with at least one real root, and that gi has no multiple roots then, 
as it is proved in appendix 0, only the coefficient Cl is non-trivial. This assumption 
turns out to be true for the cases we are studying in this paper. The boundary 
condition for the trapped surface then satisfies 



( E4 \W g,{kr) 



= 8 (6.20) 

c 



The high-energy limit of the trapped surface is defined by 

Mi » («-^^' 

(see equation (|6.21|) ). Equation (|6.20|) defines the boundary of the trapped surface 



as a function of the energy E and the wavenumber k. We will find specific examples 
where there exist r^^ = and for suitably chosen 6(r) and k in subsection |6.2.1| . 

6.2.1 Discrete spectra 

We consider a confining scale factor b which asymptotes to AdS in the UV and 
behaves as in IHQCD in the IR, 

6(r) = -e"^ (6.22) 
r 



Although the precise solution, that fits YM data is slightly different, ||3^, we will 



use the one in ( |6.22D as it does not affect the high energy asymptotics. This solution 



has been discussed in appendix G of The length scale R is the analogue of the 
QCD scale. We may parametrize it as i? ~ ^qcd- 

Unlike non-confining cases, in confining bulk geometries the physical spectrum 
of glueballs, namely the fluctuations of the metric and the scalar dilaton, is discrete, 
[EH]. The radial equations that define both the shockwave profile, in equation 



as well as the equation for the trapped surface, (|3.8|) involve the same radial equation 



as the one that determines the spectrum of spin-2 fluctuations (2++ glueballs). 

When it comes to determine the presence of a trapped surface, a similar condition 
appears: there is no trapped surface unless the solutions to ( p. 81) are normalizable 
(see ( |6.24| )). This is directly obvious in Shockwaves determined by a single graviton 



wavefunction and can be shown in general. 

As normalizability of the wavefunction is important for the existence of the 
trapped surface, the spectrum of shockwave proffies is discrete, with the transverse 
momentum determined from the 2+"*" glueball masses m„ by \kn\ = mn- Therefore, 
the transverse profiles of the shockwave distributions, not surprisingly, are deter- 
mined by the normalizable 2"*"+ (graviton) wavefunction. 

For the metric (|6.22|) , when k takes the particular subset of values 



12 

kl = {n + 2)— n = 0,1,2,3,... (6.23) 
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the radial solution of equation reduces to a (finite) polynomial that behaves as 
at small r and is normalizable. Normalizability in the radial direction is defined 
when the corresponding eigenfunction 5'i(A;„r) satisfies 



J b{r)^\gi{knr)\'^dr < oo. 



The set of the normalizable eigenf unctions is given by 



i?4 



L^fHSrV^'), ^ = 0,1,2, 



(6.24) 



(6.25) 



where Ln^ are the (finite) associated Laguerre polynomials of degree n. This is the 
qualitative behavior in any background that is confining with a discrete spectrum of 
glueballs and a mass gap. The values in ( |6.23D coincide with the mass-spectrum of 
2++ glueballs. 




R r 

Figure 1: The (closed) boundary x^c{fC^E) of the trapped surface for the first excited 
state corresponding to 6(r) = L/re~^ for fixed energy E. There is a family of such 
curves; one for each value of E which moves upwards for larger E. The section x_l = is a 
part of the boundary as should because the source has coordinate x_l = 0. The surface does 
not extend beyond r > R which implies that there is no entropy production for energies 
greater than a scale which is naturally identified with Aqcd- 



We now analyze the example ( p.22| ) and study the formation of a trapped surface 
via collision of such quantized transverse distributions. The case n = gives a shock 
that behaves as gi ~ for all r and hence it has only a single root (see discussion 



after ( |6.21[ ). We therefore conclude that this ground state mode will produce no 
trapped surface and therefore no thermal medium. 

We next considering the n = 1 case which yields the two independent solutions 

„2 \ 



g^ihr) = (36)- 



g2{kir) = 



r / r 



i?2 



(6.26) 
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The solution gi has the right asymptotics in order to generate a closed trapped surface 
for all energies. The shock given by gi has r'^ behavior at small r in accordance with 
the expectation value of the gauge theory stress-energy tensor. The boundary along 
the r direction is given by rc^ = and = R- Indeed, when r takes these two 
values, the numerator of the left-hand-side of ( |6.20| ) becomes zero and compensates 
the zero of the denominator when = as figure ^ depicts. 

The (two pieces of the transverse) exponential profile trapped surface (without 
any cut-offs) yield 

stap = -2T7u--7r, / — r ^(^' y)dy. (6.27a) 



ERKl\g^{y^l2{n + 2))\ 



x±Io{x^) (6.27b) 



87rLV24(n + 2) 

where yo is the highest root of gi{y 12{n + 2)) = and n denotes the n^'^-excitation. 
For the case at hand, n = 1 while gi is given by (|6.26| ) with y = r/R. 

It is evident that the entropy, unlike ||31|, depends not only on the transverse 



size R/ y 12(n -|- 2) but in addition on the the confinement scale 1/R. In fact, these 
two parameters may be varied independently. 

We may compare with the analogous AdS calculation where the transverse profile 
is taken to be the same as here, with characteristic scale k. The confining theory 
has another scale 1/R that can be traded with varying the integer n. We find that 
for any n the area of the trapped surface in AdS is always larger than that in the 
confining background. 

7. Accommodating asymptotic freedom 

As was pointed out in |3^, the UV (small r) should not contribute importantly to 



the entropy production S. The reason is that by definition for perturbation theory 
to be valid, the generated particle multiplicities must be small. Many examples are 
known when large multiplicities imply the breakdown of perturbation theory, with 
the sphaleron case the most prominent one, 



We expect that in the QCD UV, at some point perturbation theory takes over, 
and this is defined as the regime in which the generation of entropy is small compared 
to that generated from lower scales. This transition we will approximate as an abrupt 
transition: we will assume that this is a radial position ruy, below which we can use 
the gravitational description, while above it standard perturbation theory takes over. 
We will neglect the perturbative contribution to the entropy as we expect it to be 
small. 

Therefore in this approximation we will introduce a UV cutoff r = rjjv iii the 
trapped surface that will simulate the emergence of weak coupling in the UV. The 
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position of ruv must be determined, and at this stage it appears as an additional 
phenomenological parameter. 

A related question is to what extend geometries with varying coupling constant 
like IHQCD implement the fact that interactions are weak near the UV, as we would 
have expected from QCD. The answer is that in asymptotically AdS backgrounds, 
even as the string coupling ^ in the UV, the graviton interactions remains 
strong. It would have been probably different if the geometry becomes AdS in the 
string frame. However in this ) the boundary geometry is singular (b) Such 

a case cannot be a solution to a gravitational action with two derivatives only. 

In the sequel we will cutoff the geometry in the UV and explore the result. 

7.1 Energy-independent cut-off 

The work suggests that in the high-energy limit, the entropy, for the geometries 



studied, should be given by 

StrapiE) ~ / b{rfxl{E ^ oo,r)dr, rc,{E) -> C, (7.1) 

where Ci is an energy independent constant^^. When the energy dependence of x± 
is of the form x\ = J2i Ci(-E)c^(r), as it is usually the case, the last equation reduces 
to 

StrapiE) ^xl{E^oo,r = Ci). (7.2) 



In particular, the discussion above and equation (3/7) imply that the geometries 
( |6.1.1| ) with a UV (constant) cut-off yield 

2 1 

Strap ^ qc ^ ~ si+3a a > 1/3, (7.3a) 

Strap qc ^ ~ a < -1. (7.3b) 

We remark that for a > 1, the UV cut-off is placed on the upper bound of the integral 
in ( [7. 1|) as the boundary theory is at r = oo in this case. Evidently, this procedure, 
modifies the (center of mass) energy dependence (that we denote by s) of S. In the 
case of [^, it reduces to S" ~ se from ~ ss (see equation ( |6.7| ) and second line of 
table ^ for a = —1). 

7.2 Energy-dependent cut-off and the saturation scale 

We will now consider the shocks with uniform transverse dependence for simplicity. 



^^In [0 the authors consider in addition (E) — >■ C2 in order to remove the IR contribution 
to the entropy production. In our case, this has aheady been taken into account by the confining 
IR geometry. 
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In all examples we have analyzed, in the high-energy limit the trapped surface 
produces the following entropy 



trap 



r)dr 



s — )■ oo 



(7.4) 



rc{E) 



where bQ{r) is the asymptotic form of b for small r. rc = rc{E) is the (lower) 
boundary of the trapped surface and it is determined by the boundary conditions 
(see e.g. (|6.20|) ). It is generically energy-dependent. We believe, that this behavior 
is much more general. 

Equation (773) implies that most of the entropy originates in the UV part of the 
trapped surface. As the energy becomes larger this part enters into the region of 
the asymptotic freedom where the coupling is small and where we do not expect a 
large multiplicity to be produced. Therefore, as we have argued in the beginning of 
this section, we will impose asymptotic freedom by cutting-off the trapped surface 
at some tq as in [0. It is natural to expect that ro may be energy-dependent. 
We propose as a natural cut-off the saturation scale Qs (see figure H) by identifying 



I.e. 



1/0,, 



R 



Asym. Freedom Entropy Production 



IR 



Figure 2: The entropy production occurs between the scales Aqcd ~ ^/R and Qs which 
in the geometric language corresponds to l/Qs < r < i2. The Aqcd has akeady been 



incorporated by the usage of a confining scale factor b{r/R) as in ( 6.22 ). The weak cou- 
pling regime on the other hand can not be implemented rigorously in the framework of the 
AdS/CFT duality because string theory is not weakly coupled at high energies where 't 
Hooft coupling is small. The phenomenological Qs cut-ofi' is proposed as a first approxi- 
mation. 



5'frap/trans. area ~ / fop (r)(ir and 5 ~ / b^{r)x']_{s ^ oo,r)dr, s oo. 
Ji/Q, J HQ. 



'l/Qs Jl/Qs 

(7.5) 
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In (|7.5| ) the scale Qs can be computed perturbatively and may be modeled by 

Ql{E) ^ {Q.2GeVf x A^l'^^^f^ (7.6) 

where ^/snn is measured in GeV and denotes the cm. energy in nucleon-nucleon 
collisions. In this equation A is the atomic number while typically A lies in [0.1, 0.15] 
for energies at RHIC (and LHC) depending from the nature of the nuclear matter (pp, 
AA etc) participating in the collision |3|, ||, |5], |6|, 0, 0, |2|, ||, g 



and ^/smi is measured in GeV's. The interval for the values of A is obtained from 
fitting data from independent processes, e.g Deep Inelastic Scattering (DIS) while 
they are close to the results predicted by analytical calculations . In the uniform 



transverse dependence case and assuming that at small r the geometry maps to AdS^ 
then cutting the surface at 1/Qs yields 

Strap ^ Ql S ^ OO. (7.7) 

It is interesting to note the dependence for any scale factor 6(r) that has the 
asymptotic behavior bo = L/r\ The same dependence of entropy with was found 
recently in for uniform transverse nuclei in the context of a different approach. In 



addition, multiplicity densities dN/drj are proportional to ~ as these are predicted 
on theoretical grounds from perturbative methods r] denotes the pseudorapidity 
while dN/dr] ~ Ql describes the data quite well [|7|, H H- 

However, when non- uniform transverse dependence is present, we can not reach 
such a general conclusion as the one provided by { \l.7\ j (compare with ( p.4| ) for in- 
stance). 

8. Results 

8.1 Behavior at Large s 

Using the asymptotic form of the boundary conditions, equations ( |6.7|) , (|6.13|) , and 
equations (|3.7|) , (|6.27| ), ( [7. ID , ( [7.2|) and (|7.5|) wherever appropriate we analytically 



or numerically compute the entropy production for all the cases we have considered. 
The results are summarized in tables |], ^ and ^. 
We remark that: 

• (a) It makes sense to cut the surface at ~ l/Qs if and only if > r^ placing 
restrictions on the allowed values of a in the power-low case (see caption of table 
I and footnote 16). 

• (b) The power-law shocks for a = —1 (see [0) when the surface is cut at 
1/Q, yields Strap ~ r'Qs{r'Ay/^s^/^ ~ ^17/1830.24 j^j^. ^j^^^^ ^/ 

is the radius of the nucleus A. This energy dependence describes data well (see 
plots of figure In addition, an almost linear dependence with the number 
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of participants (see figure in PbPb collisions at 2.76 TeV has been observed 
|63| . The result for 2.76 TeV concerns the ALICE experiment at LHC and is, 



up to this time, a preliminary result (see figure 1). 

(c) The case of ( |6.22|) has been analyzed for the exponential (in the transverse 
direction) profile corresponding to ki. The results found fit satisfactorily the 
RHIC data^^ 2OO GeV (see plots of figures |). 



(d) The geometry corresponding to ( |6.10D is also interesting. In this case a 



trapped surface may be formed when a cut-off is placed in the UV region. 
The resulting entropy then behaves as Strap ~ log^(s) (see 0, H) at large s 
even when the cut-off of the trapped surface is energy- dependent. We remind 
the reader that the geometry (|6.10|) is a "marginal" case corresponding to 
continuous spectra with a mass gap, 

(e) We have numerically derived a set of inequalities about the entropy pro- 
duction. In all the cases analyzed and compared (see below), the energies are 
taken identical and large^^: 

1. For an exponential transverse profile corresponding to the geometry of 
we find ^tap > Silp > ^tap > ER = fixed. Taking into 



account that kn sets the transverse size, we conclude that more dilute 
transverse distributions at fixed energy, result in more entropy. 

For an exponential transverse profile corresponding to the confining ge- 
ometry of ( |6.22D we find that as ER — )■ 00 with E/k fixed, the entropy 
increases and becomes that of the geometry b = L/r (AdSs) (see ( pTTP ). 

For the geometry b = L/r we compare the (trapped) entropy of a shock 
with exponential (transverse) profile (exp shock) with a shock having a 
power-like transverse profile, |^ (GPY shock). We assume that both of 
the transverse profiles fall-off for x± > l/k; the first falls-off exponen- 
tially (as KQ{kx±)) while the second as a power, l/i^x^ + 1/fc^)^. We find 
^irap' > Strap ^e couclude (again) that more dilute energy distribu- 
tions produce more entropy at the same total energy. 

For uniform transverse distributions, the analysis of the geometries b = 

L/r (AdSs) and ( 1021) results in Sf/^^' > S^ap'^'^'"' ■ We conclude that 
confined matter produces less entropy than conformal matter at infinite 
coupling. This is accord with basic intuition. We point out however, that 
the difference in the two entropies is subleading at high energy because 
most of the entropy is produced at the UV where the two geometries 
coincide {L/re~^^^^^ ^ L/r for small r.). 



"'^'^See next subsection under which circumstances the fitting of (any) data is achieved. 
"'^^We do not assume a (UV) cut on the surface(s) in deriving the inequahties. 
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Table 2: Classification of trapped surfaces and entropy production S = Strap at high- 
energies of shocks with uniform transverse distribution. In aU the cases, the last column 
assumes that 6(r) reduces to L/r at small r. The first column of the table displays the 
large r asymptotics. The cut-off at the last two columns refers to a cut-off of the UV 
region of the trapped surface. 



b(r) 


Confining 


Sources 


S with 
no cut 


S with 
cut at 
const. 


S with 
cut at 

i/Qs 




For a > 


No 


~ s-*-/^ log(s)"'^^ 


~ s^/^ log(s)"'^^ 


-Ql 


r R \a 

(ro - r)" 


For a > 
For a > 1 
For a < -1 


No 
No 
No 


~ s^/^ log(s)^ 

3a+l 
~ S 6a 

3a+l 
~ S 6a 


~ S"*"/^ log(6")^~ 

3a+l 
~ S 6a 

3a+l 
~ S 6a 


-Ql 



Table 3: Classification of trapped surfaces and entropy production at high-energies of 
shocks with non-trivial transverse dependence. It is assumed that ~ with A = 0.15 
(for AA collisions). 



b(r) 


Con- 


Sources 


S with 


S with 


S with 




fining 




no cut 


cut at a 
const. 


cut at 


{t-toY a>l/3 


Yes 


Yes 


3(a+l) 
rvj g2(Za+2) 


1 

rsj 5 3a+l 




(r — ro)" a < — 1 


No 


Yes 


(3a+l) 
^ 6 a 


1 

^ g-Ml-a) 


2 + :5A(;-ia2-l) 
S fi(l-a) 




Yes 


Yes 




~ log^(s) 


~ log^(s) 



8.2 Fitting Data 

8.2.1 Relating Strap with multiplicities 

The trapped surface analysis does not give the produced entropy but it provides a 
lower bound 

Strap ^ Sprod.- (^-l) 

Moreover there are several simplifying assumptions that remain between any com- 
parison of the calculations done here and experimental data. We have spelled them 
out in the introduction and commented on how much each of them is expected to 
affect the connection with the data. In particular we have argued that for collisions 
and bulk geometries where the majority of the trapped surface area comes from 
the UV, AdS part of the geometry, we should still expect that the relative factor 
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Table 4: Classification of trapped surfaces and entropy production at high-energies of 
shocks with non-trivial transverse dependence. The entries of the first line-last three 
columns is a pure guess motivated from the results of b ~ 1/r. The second line corre- 
sponds to the normalizable ki mode. The quantities m{s) and n{s) are slow functions of s 
ranging (approximately) in the intervals [0.42,0.52] and [0.5.,1] respectively as ^/s increases 
in [20, 200] GeV. In the same energy interval, ci ranges in [300, 775]. An accurate fitting is 
found (see left plot of figure @)). 



b(r) 


Confining 


Sources 


S with 
no cut 


S with 
cut at a 
const. 


S with 
cut at 

i/Qs 


L^-r/R 
r 


Yes 


Yes 


-5^/3 log2(s) 


\og\s) 




L^-r^R^ 
= 1) 


Yes 


No 


~ (^)™Wx 
(Numerically) 


Not 
Interesting 


Not 
Interesting 



relating the final entropy release to the area of the trapped surface to be almost 
energy-independent. 

We must also quantify the relation of the total multiplicity and the produced 
entropy. The total entropy is given by the the number of charged particles Nch times 
~ 3/2 to account for the neutral particles multiplied by ~ 5, that is the entropy per 



particle [^ Hence, 

Sprod. ~ 7.5N,h. (8.2) 

We have analyzed the two cases: (a) The case of the AdS geometry with a cut-off 
at the UV (at ~ l/Qs where Qg the saturations sale), (b) The case of the confining 
IHQCD-like geometry ( |6.22| ) for n = 1 (first excitation) without any UV cut-off. 

It should be stressed here that a single parameter is used in these fits, namely 
the overall constant coefficient of the leading (large) s-dependence. 

8.2.2 Multiplicities for the AdS-Q^ setup 

This setup has an AdS metric, a GPY-like transverse profile and a UV cutoff at 
r ~ 1/Qs with Qs given by ( [7.6] ). The gravity parameters are chosen according to 
[PI|] as L^/G^ ~ 1.9. We need in addition the following relations 

1 = 0.197 GeV.fm E = = (8.3) 

2 2 ^ ^ 

where A is the atomic weight of the participating nuclei in the collision and ^/SNN 
the center of mass energy/nucleon 
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We drop the subscript NN from ^^snn from now on for simplicity. 
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Figure 3: Total multiplicities as functions of the cm. energy measured in GeV. The red 



dots are experimental data from RHIC (AuAu collisions) taken from |41] with the error 
bars included while the theoretical curves are drawn using dashed lines. The same applies 
for figure^ The left plot concerns the AdS background as in |31] with a cut of the surface 
at the UV at c/Qs where c ~ 1 and transverse size r' = vau ~ 8 fm (with A= 197; see 



( |8.4| )). The right plot is the extension of the left plot at higher values of the energy for 
PbPb collisions (that is for A= 207; see (8^)) with the points for ^/s = 2.76, 5.5 and 7 
TeV inserted. 



We will cut-off the trapped surface at = c/Qs, where c a positive constant^^. 
One then computes (at large s) 



trap 



2 X 



-2tx 



lx\{r,Er' 



where Er' = 4- 



- (gy)(Ar'v^)V3 



2 -G5 
1900 ( ^'^'^^ 



+ (r — r') 

0.483 



2\ 3 



\GeV 



rAu~Sfm 



where we normalized the formula of Strap with A^u = 197 for AuAu collisions. In 
normalizing, we used the fact that r' is the transverse size of the colliding glueball^^ 
(beam) assuming that it satisfies the empirical law (r'/r^u) = A}!'^ j A^j^ which applies 
for nuclei. In the last equality, we have used ( [7.6|) with A = 0.15 for AA collisions. 



D and LVG5 = 1.9 |3T[ 



The limits of the r integration are found by setting = (see figure |T]). 
The excellent fitting of plot |] with RHIC data is achieved for K=Kau = 197 and 
( 1 /c) X (overall coefficient ^^)~ 1.54 taking tau = 8 fm^^. It is pleasing that both, the 
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As mentioned earlier, it makes sense to cut at c/Qs iff c/Qs > '"Ci = {Ay^G^r' /L^ 



-1/3 



(r' 



is the transverse nuclear size). For c ~ 1, A any (reasonable) value (see ([7.6D), r' « (A/A^u)^/"^ x 8 



fm and > 20 measured in GeV, this condition is satisfied. 

^^It corresponds to the position of the point-like source in the fifth dimension (see discussion 



above (|6.5|)). 

^*The meaning of the overall coefficient is discussed in subsection B.2.1 
^^A different value for but close to 8 fm may still fit RHIC data. 
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overall coefficient and c can simultaneously be of order one. The extrapolation to 
higher energies is done in section ^]3| and the right plot of figure ^. 

8.2.3 Multiplicities for the IHQCD setup 

This setup involves a metric b = L/re~^^^^^ without a UV cut-off and an exponential 
transverse profile. In this case, Strap is given by ( |6.27|) for n = 1 and has to be solved 
numerically. The parameters of gravity were determined by matching lattice data in 
[p9|] and we will use these in the sequel. Choosing A'^c = 3 we find 



6 

ki = - = 3.1 GeV — ^ 1.96, ( 

R ACg 

where ki = mi is the mass of the lightest spin- two glueball (second excitation). 



.5) 
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Figure 4: Both figures concern the numerical analysis of the surface (resulted from the 
first excitation) of ( |6.27| ) for lower and higher energies (dashed plots). There is not a UV 
cut-off in this case. In the left figure, the green plot is given by and exhibits the 

accurate approximation of the numerical plot found for A=A^„. The overall coefficient of 
the numerical plot has been chosen in order to fit the RHIC data. The right plot is the 
extension of the approximating plot at higher energies for A= Api^ and hence, according 
to it is given by N^h = 79.8 {^/sf■^^^ logO-'^^^ (562^5). The points for ^ = 2.76, 5.5 

and 7 TeV are inserted. 

The numerical result for multiplicities are plotted in figures ^ for weight A= Aau 
and compared with RHIC data. The left dashed plot is the result of our numerical 
analysis for A=A^„ and energies up to 250 GeV including the RHIC data. The 
agreement is satisfactory. The numerical result can be approximated very accurately 
(see green plot of figure H) by 

where a/s is measured in GeV. In order to go higher in the energies, a more refined 
numerical analysis is needed. Hence, for the present work, we use the fitted curve 



- 24 - 



given by equation (|8.6|) in order to predict multiplicities for higher energies. We do 



this in section ^]3| for the energies to be reached by LHC taking into account that 
for Pb A= 207. The corresponding plot is the one on the right of figure ^. 

We remark that cutting the surface at some UV cut-off does not improve the 
fitting. In particular, for energies higher than 200 GeV we either find very low multi- 
plicities (for a constant cut-off) or very large multiplicities (for an energy dependent 
cut-off). 

We close this section by noticing that in order to obtain a (more realistic) value 
for in the case of the trapped surface given by ( p.27| ), the black hole ansatz for 



the scale factor of ( |6.22|) should be solved. We postpone a detailed analysis with the 
precise IHQCD profiles for later work and choose the particular values mentioned 
below equation (^.2|) in order to fit the data. 



8.3 Multiplicities at LHC energies 

In the previous section we have analyzed (a) the geometry of AdSs with a UV cut-off 
and (b) the first excited state of ( |6.22| ) without any UV cut-off. We have seen that 
both fit the RHIC data up to 200 GeV satisfactorily (see plots of figures ^ and ^). 
For the aforementioned geometries (a) and (b), we may extrapolate them at higher 
energy in order to assess what they predict for the multiplicities at the energies 
reached by LHC. 

• Geometry (a): Multiplying formula ( ^.4|) by 1.54 (ignoring c; see subsection 
g.2.2|) , dividing over 7.5 (see subsection ^.2.1|) and taking A = 207 for Pb 



central coUisions we find Nl^fl> ^ 18750 for^o 2.76 TeV, A^",^^ ^ 261800 for 5.5 



TeV and N^^^ ^ 29400 for 7 TeV. 

r21 



For high-multiplicity proton-proton (pp) central collisions where A= 1 we 
find N^f^ ^ 70 for 0.9 TeV, N^,^ ^ 110 for 2.36 TeV, N^^^ ^ 190 for 7 TeV 
and N^i^ ^ 260 for 14 TeV. It is pointed out that extracting experimental 
results (for total N^fJ from ATLAS Q and CMS |l], H (for these energies) 



is not a trivial task. It involves model dependent procedures and Monte-Carlo 
simulations. 



Geometry (b) The numerical analysis corresponding to ( ^.6] ) (see right plot of 
I) for PbPb cohsions predicts N^h ~ 19100 for 2.76 TeV, N^h ~ 27000 for 5.5 
TeV and N^h ~ 30500 for 7 TeV^^^ 



^°A11 the energies that are mentioned in this subsection refer either to y^SNN or to ^JSpp. 

^^The notion of centrahty in pp coUisions is trickier to define. The best definition is to select 
high muhiphcities in the final state. The definition of high multiplicity is ambiguous, but may still 
defined. 

^^A prediction for pp collisions using ( ^.6[) is less reliable for the moment because the errors 
induced from factors of Ap/AAu — 1/197 are larger than in PbPb collisions. To deal with this, a 
more refined numerical work is needed. 
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In the above results, we have assumed exactly a zero impact parameter and hence 
these are zero centrality processes. This implies that we might predict slightly larger 
multiplicities than the upcoming data. The difference between the two predictions (at 
a given value of \/s) for the two cases considered above is ~ 6% for PbPb collisions. 
It is interesting that the first preliminary results from ALICE at LHC for ^/s = 2.76 
TeV and PbPb collisions give N^h ~ 17000 ± 1000 for - 5% centralities In 
order to make a better estimation, we consider the right plot of figure 7 in p3| which 



shows that the maximum number of participants for the particular measurements, 
is approximately 380. Hence, plugging A= 380/2 = 190 and ^/s = 2.76 TeV in (|8l^ ) 
results Nch = 17300 which is extremely close to the (preliminary) measurements (see 
figure 

9. Outlook and discussion 

The initial motivation of this work was to estimate the produced entropy of colliding 
shocks in different circumstances and obtain a qualitative answer for the results, in 
view of applications to heavy-ion collisions. In particular, the goal was to compute 
Strap and hence (estimate) the particle multiplicities (see section |8.2.1| ). We have 



analyzed various dilaton-gravity geometries, with or without transverse dependence, 
with or without confinement and with or without UV cut-offs. The cut-offs were 
assumed energy dependent or energy independent. We conclude the following: 

• Both our analytical calculations and our numerical analysis have shown that 
most of the entropy comes from the UV part of the geometry provided that 
the geometry is asymptotically AdS. 

• There are geometries with a UV energy-independent cut-off that lead to an 
asymptotic ~ log^(s) behavior for the entropy (see table H). 

• For uniform transverse distributions, the AdSs geometry, produces the least 
entropy among the geometries of section This result assumes that the non- 
conformal geometries in |^ survive until the boundary. They are therefore not 
asymptotically AdS. A similar result would be valid if the associated geometries 
are asymptotically AdS, but the energies such that the trapped surfaces do not 
penetrate the AdS region. These results are summarized in table 

• We have constructed exact shock solutions, with non-trivial transverse distri- 
butions and point-like bulk sources, and have computed the entropy of the 
trapped surface (see table |). Inserting an energy-independent UV cut-off (see 
fifth column of table |^), and choosing a 1 in the first line results in 
S ~ {y/sY^'^ (see also (|7.3a| )). The particular energy dependence seems to 
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This case correspond to a power-like scale factor b(r) 
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describe RHIC data as equations (|8.4| ), ( p.6|) and figures || and ^ suggest. A 



similar power-law dependence is obtained assuming a Landau liydrodynamical 
behavior after the collisions . A characteristic of the geometries corre- 
sponding to a = 1 and —1/3 is that they do not reduce to AdSs at the UV 
unlike the a = —1 geometry. 

For geometries that have a mass gap and confinement, the entropy production, 
independent of transverse distribution, is subleading to that in AdS at the same 
total energy and transverse scale. Here it is assumed that all geometries are 
asymptotically AdS. 

The entropy production decreases when higher glueballs collide (see p.2.i|) . This 
implies that more dilute transverse energy distributions produce more entropy 
at fixed total energy {ER = fixed, where R can be identified with 1/Aqcd)- 

We will denote the trapped surface areas as follows: S^Zp^'^'^^''' for the GYP- 
like Shockwave in ( |6.5| ) in AdSs where k is the scale of the transverse profile, 
^trap^ for an exponential transverse profile in AdS5 as in ( |6.16D , and for 
an exponential transverse profile with scale /c„ for the confining metric (|6.22|) . 
In all of the above the total energy, and the transverse scale are kept the same. 

We find the following inequalities for large E: 



qGYP;AdS ^ qk;Ads „ rifcs ^ ciA;2 > ofei 
'^trap '^trap ^ ^trap r~-i ^trap ~ ^trap 



(9.1) 

fc=fel=fc2=fc3 



We conclude that the entropy production increases as the confinement scale 
(~ 1/i? see subsection |6.2.1|) decreases provided that the transverse size (with 



respect to the energy) is kept fixed. 

Equation ( p.4|) suggests that multiplicities should grow with A as Al-"^!^^ . In 
fact, an almost linear dependence of Nch with the number of participants was 
observed in the recent ALICE experiments performed at 2.76 TeV for 
PbPb collisions. According to figure H our result shows a similar behavior. 
The agreement becomes better as the number of participants increases, that is 
as the collision becomes more central which is the case that we assumed in this 
paper. We remark that the ALICE results are still preliminary. 



We have derived two formulae for multiplicities, for the AdS — Qs setup in ( |8.4| ) 
and the IHQCD setup in ( p.6[ ). We normalize the multiplicities using RHIC 
data and we then compare with known LHC data. We find that in both cases 
they are very close to the 2.76 TeV heavy-ion LHC data, |^ (see subsection 



An important puzzle of our analysis is that the scattering of the lowest lying 
2++ glueballs does not seem to lead to a trapped surface. 
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20000 



Figure 5: The dashed plot refers to our theoretical prediction for PbPb collisions at 
^/s = 2.76 TeV. It shows the total (charged) multiplicity (Nch) as a function of A (see 
(8^)) at fixed energy s. The red points are data taken from reference [^3| with the error 



bars included: the horizontal coordinate of these points is taken to be Npart/'^ where Npart 
is the number of nucleons (of the two Pb nuclei) participating in the collision at the fixed 
value of = 2.76 TeV. The agreement of our theoretical prediction with the data seems 
to improve as Npavt increases, that is as the collision becomes more central. 



For future work, we propose to use the shock of equation ( p4D and estimate the 
produced entropy. This shock has the advantage of localizing the bulk sources and 
in addition, it retains the ingredients of a mass gap and confinement. Also, more 
accurate calculations that will specify the gravity dimensionless parameter /G^ 



see e.g. (|8.5| )) could be made. Finally, it is important to understand the kind of 



black holes which give rise to the shocks (after boosting) and which in turn, mimic 
the nuclear matter in the colliders. This matter is best described by finite size black 
holes In particular, in a description of the collisions in term of what is 

known about such black holes was described. It would be interesting to make this 
picture precise, and in this numerical work will be required. 
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A. Introduction to the theory of trapped surfaces. 

The starting point is the shock-wave metric 

ds'^ = 6(r)2 {~2dx^dx- + dx\ + dr^ + 0i(a;\ x^ r)5{x^){dx^f) . (Al) 

Associated with this shock-wave we parametrize the (half piece of the) trapped sur- 
face 5*1 by 

x+ = X- + i^i(a;\ r) = (A2) 

where ipi will have to be determined. Before finding the differential equation satisfied 
by ipi it seems necessary to perform the following coordinate transformation which 
eliminates the 5{x^) from ( |A1| ). In particular we use the transformation 

x~ x~ + r)6'(x+) (A3) 

which transforms (^l|) into 

ds^ = h{rf [ -2dx+dx~ + dx\ + rfr^ - d{x^) ^ (V,0irfx^)rfx+ j . (A4) 

V j=l,2,r / 

It is (sometimes) useful to exchange the functions 0i and ipi with $i and \l'i defined 
via 

$1 = 6(r)0i *i = h{r)i)i. (A5) 

The trapped surface associated with the shock 0i can be found solving the equation 

e = h>"'V/,-^^^ = (A6) 

where a few explanations are in order: 

(a) The vector ll^^^'^^ can be generally expressed by 

l^+^Ux" = Adx+ + B{dx~ + ^#i) = Adx+ + B{dx- + 1{Y1 ^ii^i)d^') (A7) 

i=l,2,r 

which implies that 

f;'^ = e^)) = {a, B, IbV.^P,, IbV.^P,, l^V.^,) (A8) 

for some functions A and B which are determined by the following requirements: (i) 
g'"'ll:t^hi~^^^ = where gt"^ is the (inverse) metric of (^, (ii) /(+^)* > and (iii) 
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The necessity of the presence of the superscript + is explained in what foUows. 
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/(+!)- < 0. The vector l^^^^ as defined is transverse to 5*1, null, forward and outward. 
There is another vector, the which is also transverse to Si, null and forward but 
inward and can be taken to be 

= {li-'\t'\l[-'\li-'\lir'^) = (1,0,0,0,0). (A9) 

(b) The covariant derivative is taken with respect to the metric 

(c) The tensor /i^j^ is the projective metric to ll^^^ and is given by 



7(+l)/(-l) , 7(+l)/(-l) 

t/^ l-U T I'll 



(AlO) 



where g^'^ is the inverse metric of equation ( |A4D . This tensor clearly has the property 

h^^•^l{+^) = 0. (All) 

We are now in position to follow the explanations (a), (b) and (c) mentioned above. 
Explanation (a): From g^'^l^f^hl^^ = and the metric ( [A4D one finds that 



B A 



8 ^ 



^1 $1 



'6(r) 6(r)' 



6(r) h{r) 



0. 



(A12) 



i=l,2,r >- 

We emphasize that (i) the solution, 5 = does not satisfy the rest requirements of 
I'l^^'' and hence it is discarded, (ii) Both of the subscripts (j) in the Vj's of ( [A12| ) are 
lower indices. We point out that condition ( |A12| ) does not specify A and B uniquely 
but it only specifies their ratio 



A _l ^ 



i=l,2,r >- 

Hence a choice for A and B is 



'6(r) 6(r)' 



A 



i=l,2,r >- 



'6(r) h{r) 



X Vi 



^1 



h{r) h{r) 



'h{r) h{r) 
$1 



B = -i 



(A13) 



(A14) 



which is proportional (by the factor — )■ -p-) with (101) of 0] . One may 
wonder that this arbitrariness in choosing A and B may invalidate the procedure 



that determines the trapped surface and in particular equation ([A6|) . We show later 
that this is not the case and that ( [A6| ) is invariant under rescalings of the vector 
Now, according to ( [A8| ) , the choice ( [A14| ) for A and B implies 



/(+!) _ f/C+i) /(+!) /(+!) /(+!) /(+!) 

^1 




V,;( 



$1 



6(r) h{r) h{r) h{r) 



,4Vi( 



^1 
6(r) 



,4V2(7T^),4V/ ' 



'6(r) 



h{r) 
(A15) 



^^Up to a factor of 2 for B which is traced to the authors different convention for the hght-cone 
coordinates. 
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Explanation (b): One has then to compute the covariant derivative of /^^^^ with 



respect to the metric (|A4|) forming the tensor V;//^^^''. We do not display explicitly 



this tensor. 

Explanation (c): The final step is to contract 'V,yll^^^ with /i^"^. This yields 

(□^,53-mr)))(^i-<fi) = Wr)) = ^^^^^1^^ (A16) 

where nAdSa defined with respect to the metric 

ds^ = h{rf [dxl + dr^) . (A17) 

Associated with a second shock 02 moving along x~ there exists the surface 
5*2 constructed from a second vector the analogue of Z^^^'*. From symmetry 

considerations, I'lt'^'^ can be immediately found from /|t^^^ (see ( A15|) ) and hence 
reads 



/(+2) 




V,: 



^5 



$9 



h{r) h{r) h{r) h{r) 



,8,4Vi( 



^2 

h{r) 



, ^2 



),4V, 



._*2_ 

'6(r) 



(A18) 



Boundary conditions on C = fl 6*2: 
\E'i|c = ^1^2 1 c (continuity of S on C) and ^ 



i=l,2,r 



E 



i=l,2,r 



• • ;(+l) 

requirmg 



and Vi^'i 



C — 
l,2,r (by 



c 



I 



(+2) 



$2 = $ and by symmetry \E'i 
we get finally are 

1 



c). Assuming for simplicity identical shocks, we have $1 = 
\E'2 = ^ = on C. Hence the boundary conditions 



^1 



c 







i=l,2,r 



h{r\ 



-Vi^ - V,; 



-V,:^ - V,: 



6(r 



(A19) 



Equation ([A19| ) may be simplified. Taking also into account that in the pre-coUision 
time {x^ < 0) $ vanishes (see ( |A^ ) for s+ < 0), equation ( [A 191 ) reduces to 

1 







6(r) 



5^ [V.W.^]|c = 8. 



(A20) 



i=l,2,r 



Finally, we prove that (|A6D is invariant under rescalings of /^^^^ by an arbitrary 



(differentiable) function a{x''): We consider equation ( |A6|) with /li^^-* — )■ 0(0:'^)/)^ 



7(+l) 



We have 6 



(+1) 



(+1)^ 



Using the Leibniz rule of the 



covariant derivative we get 6' = ah^'^V^{lu ) + hf^^ll d^{a). Using the projective 



property of h'^" on equation ( |A11|) , we find that the second term of last equality 



vanishes yielding 6' = a/i^^V^(/ 



(+ih 



a6 = since = by assumption □ 
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B. Localizing the bulk sources 

The shocks we have considered in subsection correspond to a non-locahzed bulk 
source in the r direction (see (|6.3|) ). In order to locahze this source we begin from 
( |6.3|) and consider appropriate hnear combination of solutions. This is not hard as 
one merely has to use the completeness relation of Qk which schematically has the 
form 



dkC{k)gk{r)gkir') = 6ir-r'). (Bl) 
This would imply that the shock 0^ should be given by 

dkC{k)Ko{kx^)gk{r)gk{r'). (B2) 
In particular the completeness relation 

y — — -x'"e-"'Li") (x (x) = 6(x - x') (B3) 

^ r(n + a + 1) n \ J n J\ J V J V / 

n=0 ^ ' 

for the Laguerre polynomials implies that the correctly normalized shock^^ should 
be given by^^ 



00 

X 



5 (t v/^^) 4^' (4) (4) ^ (B^) 



In arriving to ( P4|) we have used 



n=0 



n\2 



C. Proof of formula (^TZUj) 

The idea is to note that the boundary of the surface should contain the point 
{rc,x±c = 0) because the source is located at x±c = while both the terms {dripk) 
and {dx^'>pk) should be finite at x^c- We also state the following relations 

{d,,Ko{kx^))h{kx^) - K^{kx^){d^Jo{kx^)) = -— + {k' - k)Ok'{x^), (Cla) 

^^Such that the gauge tensor T++ integrates to E. 

^^For the AdSs geometry, the sources are locahzed choosing C(fc) — kJ2{kr'). The shock then 
becomes that of This can be shown using resuhs of reference pa]. 
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lim Io{kx±] 



l + 0{x] 



Beginning from {dx^^ipk) we have that 



hm d^^ipk 



r 



x±c 



(Clb) 



(C2) 



where we have used ( |(J1|) in order to simphfy the last expression. It is evident that 
9i{krc) should have (at least) two distinct real roots in order to have a trapped 
surface (see figure |I|). For reasonable shocks, the gi decays for small r as and 
hence the one root is at tq = 0. By assumption of the claim (see subsection 
there exists one more real root call it r^j and as a result the following is true 



lim gi{krc) = 



(C3) 



We now compute dripk again for x±c and hence at Tc — >■ rc^ ■ We have 

Ko{kx±) 



lim dri^pk 



L3p 



lim 



X 



E loik'x^) [g[{krc) {Cl,g,{k'rc) + ClMk'rc)) 



- g,{krc) {Clg[{k'rc) + Clg'^ik'rc)) ] 
Koikx±)g[{krc) 



E4 



lim 



(C4) 



where in the second equality we have used (p3|). The simplified expression for dr{ipk) 
implies a logarithmic divergence when x± because limr^-^rc^ 9i{ki^c) 7^ by the 
hypothesis of the claim that there are no multiple roots for gi'^^. The fact that the 
right-hand side of ( p^) is independent on the arbitrary coefficients Cl, C| would 
yield to the naive conclusion that the whole second term of ( 5.19|) should be absent. 
However, looking more carefully inside the sum of the numerator of ( p^) , one realizes 
that the terms for gi{kx±c) = gi{k'x±c) do not actually participate in the sum and 
hence they are allowed. Thus the trapped surface equation is given by 



ipk = 4>k{r,x±) - 



gi{kr)Io{kxj_) 
9iikrc)Ioikx±c) 



4>k{rc,x±c) 



yielding to the trapped surface boundary determined by the condition (|6.20|) com 
pleting the proof of the claim. 
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With the exception for r = which is a fourth- root. 
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